Probabilistic fatigue methodology for six nines reliability

Fleet readiness and flight safety strongly depend on the degree of reliability that can be designed into rotorcraft flight critical components. The current U.S. Army fatigue life specification for new rotorcraft is the so-called six nines reliability, or a probability of failure of one in a million. The progress of a round robin which was established by the American Helicopter Society (AHS) Subcommittee for Fatigue and Damage Tolerance is reviewed to investigate reliability-based fatigue methodology. The participants in this cooperative effort are in the U.S. Army Aviation Systems Command (AVSCOM) and the rotorcraft industry. One phase of the joint activity examined fatigue reliability under uniquely defined conditions for which only one answer was correct. The other phases were set up to learn how the different industry methods in defining fatigue strength affected the mean fatigue life and reliability calculations. Hence, constant amplitude and spectrum fatigue test data were provided so that each participant could perform their standard fatigue life analysis. As a result of this round robin, the probabilistic logic which includes both fatigue strength and spectrum loading variability in developing a consistant reliability analysis was established. In this first study, the reliability analysis was limited to the linear cumulative damage approach. However, it is expected that superior fatigue life prediction methods will ultimately be developed through this open AHS forum. To that end, these preliminary results were useful in identifying some topics for additional study

Multiphoton Absorption of Myoglobin–Nitric Oxide Complex: Relaxation by D-NEMD of a Stationary State

ABSTRACT: The photodissociation and geminate recombination of nitric oxide in myoglobin, under continuous illumination, is modeled computationally. The relaxation of the photon energy into the protein matrix is also considered in a single simulation scheme that mimics a complete experimental setup. The dynamic approach to non-equilibrium molecular dynamics is used, starting from a steady state, to compute its relaxation to equilibrium. Simulations are conducted for the native form of sperm whale myoglobin and for two other mutants, V68W and L29F, illustrating a fair diversity of spatial and temporal geminate recombination processes. Energy flow to the heme and immediate protein environment provide hints to allostery. In particular, a pathway of energy flow between the heme and the FG loop is illustrated. Although the simulations were conducted for myoglobin only, the thermal fluctuations of the FG corner are in agreement with the large structural shifts of FG during the allosteric transition of tetrameric hemoglobin

Fractal Analysis of Protein Potential Energy Landscapes

The fractal properties of the total potential energy V as a function of time t are studied for a number of systems, including realistic models of proteins (PPT, BPTI and myoglobin). The fractal dimension of V(t), characterized by the exponent \gamma, is almost independent of temperature and increases with time, more slowly the larger the protein. Perhaps the most striking observation of this study is the apparent universality of the fractal dimension, which depends only weakly on the type of molecular system. We explain this behavior by assuming that fractality is caused by a self-generated dynamical noise, a consequence of intermode coupling due to anharmonicity. Global topological features of the potential energy landscape are found to have little effect on the observed fractal behavior.Comment: 17 pages, single spaced, including 12 figure

Ligand binding and conformation change in the dimeric hemoglobin of the clam Scapharca inaequivalvis.

The reaction with carbon monoxide of the cooperative dimeric hemoglobin from Scapharca inaequivalvis has been examined by flash photolysis. In the nanosecond time range, geminate rebinding of 5% of dissociated CO occurs with a rate constant of 1.4 x 10(7) s-1. There is a change in absorbance of deoxyhemoglobin following photolysis at a rate of 1.2 x 10(6) s-1, consistent with a shift in the position of the Soret band to longer wavelengths. The amplitude of the change is proportional to the population of deoxydimer. In much of the Soret region this change is greater than the absorbance excursion associated with geminate recombination. There is at least one other slower change associated with the singly liganded species. Geminate rebinding of NO has components of 50, 8, and 0.035 ns-1, accounting for 75%, 25%, and less than 1% of the total reaction observed after a 35-ps photolysis flash. Simulation of diffusion of NO by molecular dynamics shows the ligands moving from the heme pocket to a subsidiary space between the edge of the heme and the surface of the protein

Action-derived molecular dynamics in the study of rare events

We present a practical method to generate classical trajectories with fixed initial and final boundary conditions. Our method is based on the minimization of a suitably defined discretized action. The method finds its most natural application in the study of rare events. Its capabilities are illustrated by non-trivial examples. The algorithm lends itself to straightforward parallelization, and when combined with molecular dynamics (MD) it promises to offer a powerful tool for the study of chemical reactions.Comment: 7 Pages, 4 Figures (3 in color), submitted to Phys. Rev. Let

Constant amplitude and post-overload fatigue crack growth behavior in PM aluminum alloy AA 8009

A recently developed, rapidly solidified, powder metallurgy, dispersion strengthened aluminum alloy, AA 8009, was fatigue tested at room temperature in lab air. Constant amplitude/constant delta kappa and single spike overload conditions were examined. High fatigue crack growth rates and low crack closure levels compared to typical ingot metallurgy aluminum alloys were observed. It was proposed that minimal crack roughness, crack path deflection, and limited slip reversibility, resulting from ultra-fine microstructure, were responsible for the relatively poor da/dN-delta kappa performance of AA 8009 as compared to that of typical IM aluminum alloys

A Doubly Nudged Elastic Band Method for Finding Transition States

A modification of the nudged elastic band (NEB) method is presented that enables stable optimisations to be run using both the limited-memory quasi-Newton (L-BFGS) and slow-response quenched velocity Verlet (SQVV) minimisers. The performance of this new `doubly nudged' DNEB method is analysed in conjunction with both minimisers and compared with previous NEB formulations. We find that the fastest DNEB approach (DNEB/L-BFGS) can be quicker by up to two orders of magnitude. Applications to permutational rearrangements of the seven-atom Lennard-Jones cluster (LJ7) and highly cooperative rearrangements of LJ38 and LJ75 are presented. We also outline an updated algorithm for constructing complicated multi-step pathways using successive DNEB runs.Comment: 13 pages, 8 figures, 2 table

Efficient Dynamic Importance Sampling of Rare Events in One Dimension

Exploiting stochastic path integral theory, we obtain \emph{by simulation} substantial gains in efficiency for the computation of reaction rates in one-dimensional, bistable, overdamped stochastic systems. Using a well-defined measure of efficiency, we compare implementations of ``Dynamic Importance Sampling'' (DIMS) methods to unbiased simulation. The best DIMS algorithms are shown to increase efficiency by factors of approximately 20 for a $5 k_B T$ barrier height and 300 for $9 k_B T$, compared to unbiased simulation. The gains result from close emulation of natural (unbiased), instanton-like crossing events with artificially decreased waiting times between events that are corrected for in rate calculations. The artificial crossing events are generated using the closed-form solution to the most probable crossing event described by the Onsager-Machlup action. While the best biasing methods require the second derivative of the potential (resulting from the ``Jacobian'' term in the action, which is discussed at length), algorithms employing solely the first derivative do nearly as well. We discuss the importance of one-dimensional models to larger systems, and suggest extensions to higher-dimensional systems.Comment: version to be published in Phys. Rev.

Two-dimensional visibility charts for continuous curves

This paper considers computation of visibility for two-dimensional shapes whose boundaries are C1 continuous curves. We assume we are given a one-parameter family of candidate viewpoints, which may be interior or exterior to the object, and at finite or infinite locations. We consider how to compute whether the whole boundary of the shape is visible from some finite set of viewpoints taken from this family, and if so, how to compute a minimal set of such viewpoints. The viewpoint families we handle include (i) the set of viewing directions from infinity, (ii) viewpoints on a circle located outside the object (for inspection from a turntable), and (iii) viewpoints located on the walls of the shape itself. We compute a structure called a visibility chart, which simultaneously encodes the visible part of the shape's boundary from every view in the family. Using such a visibility chart, finding a minimal set of viewpoints reduces to the set-covering problem over the reals. Practical algorithms are obtained by a discrete sampling of the visibility chart. For exterior visibility problems, a reasonable approach is to compute an almost-optimal solution (in terms of number of viewpoints), which can be done in almost-linear time. For interior visibility problems, or when a more correct solution is required, we solve the general set-covering problem, guaranteeing an optimal solution but taking exponential time

Constraint methods for determining pathways and free energy of activated processes

Activated processes from chemical reactions up to conformational transitions of large biomolecules are hampered by barriers which are overcome only by the input of some free energy of activation. Hence, the characteristic and rate-determining barrier regions are not sufficiently sampled by usual simulation techniques. Constraints on a reaction coordinate r have turned out to be a suitable means to explore difficult pathways without changing potential function, energy or temperature. For a dense sequence of values of r, the corresponding sequence of simulations provides a pathway for the process. As only one coordinate among thousands is fixed during each simulation, the pathway essentially reflects the system's internal dynamics. From mean forces the free energy profile can be calculated to obtain reaction rates and insight in the reaction mechanism. In the last decade, theoretical tools and computing capacity have been developed to a degree where simulations give impressive qualitative insight in the processes at quantitative agreement with experiments. Here, we give an introduction to reaction pathways and coordinates, and develop the theory of free energy as the potential of mean force. We clarify the connection between mean force and constraint force which is the central quantity evaluated, and discuss the mass metric tensor correction. Well-behaved coordinates without tensor correction are considered. We discuss the theoretical background and practical implementation on the example of the reaction coordinate of targeted molecular dynamics simulation. Finally, we compare applications of constraint methods and other techniques developed for the same purpose, and discuss the limits of the approach